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ABSTRACT 



The buckling action of hollow cylinders, subjected to internal 
pressure, was investigated under various conditions of axial load in the 
cylinder walls# The theory of Suler was found to apply where the modulus 
of elasticity and second moment of oross-seotional area are determined 
solely by the material and cross-section of the cylinder, and the load 
is the total axial load on the oolumn# 
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SYMBOLS 



Ac 


Cross-sectional area of the cylinder (in.^). 

2 

Cross-sectional area of the piston (in. ). 


a 


Lateral deflection of the mid-point of the column 
( #001 in. ) . 


L 

r 

l 

-X 


Modulus of elasticity of the material of the 
cylinder (psi). 

Second moment of cross-sectional area about a 
centroidal axis (in. 4). 


J? 

M 

P 


Length of the column (in.). 
Bending moment (in.-lbs.). 
Axial Compressive load (lbs.). 
Critical load (lbs.). 


A 

(U 


Internal pressure (psi.). 

Critical internal pressure (psi.). 
Radius of curvature (in.). 


e 


Radius of gyration (in.). 
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I Introduction 



1. Introduction, 

A lon£ cylindrical shape subjected to a compressive load will fail 
n suddon bon 1 in , provided the member is slender. This is the type of 
load.n^ invest 'ated by ,uler [" 1 } who determined a "critical load" such 
V <\t an load ’•reat^r than this will cause the column to bond. In his 
anal sis 'uler oonsidored the column to be solid and subjected to an axial 
compressive load applied at the ends. 




i -ure 1, Schematic doprosert at ion of Column 



Arsine a hollow c /Under filled with a flaid capable of boin^ 
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the length of the column in the hinged-end oonf igurat ion ; then 




A free body diagram of the column of Fig, 1 is identical to the free body 
diagram of any column with the same external restraints. Applying Eulers 
equation to the column of Fig, 1 results in: The load being the product 

of the piston areay^ p and the internal pressure^* The modulus of 
elasticity and the seoond moment of the oross-sect ional area are deter- 
mined solely by the material and the cross-section of the cylinder, 
Writing Euler’s equation in terns of a hollow cylinder filled with fluid 
results in 




In the undefleoted condition the entire axial load of the oolumn pf 
Fig, 1 is carried by the fluid. Thus we are led to the interesting 
paradox of a oolumn, which in the undefleoted condition sustains no axial 
load, buckling when the internal pressure reaches a certain oritioal value. 
Interest was focused on this problem when a long cylinder used in the 
hydraulic system of a submarine was observed to buckle when the piston 
approached the end of its stroke. With the increase in operating pressure 
of hydraulio plants and the associated decrease in size of ocmponent 
parts, it is oonoeived that in some applications the buokling aotion of 
cylinders subjected to internal pressures will be important, 

2, Object and Soope 

The object of this investigation was to verify the applicability of 
the theory of Euler to hollow cylinders subjected to internal pressure. 

It was previously shown that 
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Apt 



This equation states that, for a given length, the internal pressure which 
will cause columns of the same material to buckle is proportional to the 
quotient pf the second moment of cross-socsional area by the area of the 
piston. 

The axial load carried by the cylinder is determined by the relation 
between the area of the piston and the area of the inside of the cylinder. 
The area of the inside of the cylinder does not appear in the equation. 
Therefore it may be stated that, theoretically, the critical pressure is 
independent of the "load-ratio • n (Load-ratio is defined as the ratio of 
the axial load carried by the cylinder to the total axial load on the 
column. ) 

To verify this conclusion it was planned to test three sets of 
columns. In the first set the ^.nner area of the cylinder was less than 
that of the piston (a positive load-ratio). In the second set the areas 
were equal (a zero load-ratio). In the third set the inner area of the 
cylinder was greater than that of the piston (a negative load-ratio). 

To check the validity of the theory through a range of pressures, 
the lengths of the column were changed so that slenderness-ratios from 
100 to 200 were obtained. The effect of these variables was studied 
experimentally and the results correlated with theory. 

3 0 Acknowledgments. 
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II Theory 



Consider a straight cylindrical tube of homogeneous material, which 
obeys Hooke’s Law, filled with a fluid capable of being placed under 
pressure. One end of the tube is rigidly sealed and the other end is 
sealed by means of a frictionloss piston. The two ends are subjected to 
an externally applied axial load# The ends are restrained against lateral 
movement but have rotational freedom# The free body diagram of the column 
(Fig. 1) is identical in form to the one considered in Euler’s classic 
work# 

From Euler’s solution it is known that under an axial compressive 
loading of magnitude P, applied at the centroids of the end cross-sections, 
the tube will remain straight until P reaches a certain critical value# 

As P exceeds this critical value the tube will begin to bend laterally in 
a plane containing the axis of the tube. If the tube has a circular cross- 
section this plane will not be known in advance# The lateral bending 
increases rapidly with P, so that it is practical to consider that the 
critical value of P represents the maximum load that the column can 
support. In order to evaluate this critical load we choose a set of 
cartesian coordinates, with the x-axis the locus of centroids of the cross- 
sections of the undeformod tube, the xy-plane containing the axis of the 
deflected tube, and the z-axis perpendicular to the x and y axes# The 
origin is at one end of the tube# If we consider the column in a deflected 
position (Fig. 2) it is clear that the magnitude of the bending moment at 
a particular cross-section is Py where the y-ooordinate is that of the 
centroid of the section# 
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The condition of particular interest > the one in which the column first 
starts to deflect • At this time the slope is snail and ne^lectin^ it 
will produce a ne^li^ible error and greatly 



simplify the solution 



Thus we have 



_L 



j y 



Employing a consistent sot of algebraic signs we are led to the differen- 
tial equation 



+ 



p 



d* 



El 



y = 



where 
^ " £ Z 



The general solution of this equation is 

Z b X B fax y 

The end conditions are satisfied if 

A - 0 

and either: B --0 

or l\J*0 

If B'O then y=0 and this is the solution prior to buoklingi 
If kf-0 then 



or 



or 



kJl s /n TT, 

*(£f r 



r ^ V* £1 where /n is an integer# 



The solution of principal interest is the one obtained for , 

from which 



B * 



tt* Cl 



(The solution obtained for sn 9 0 is that of the undeformed column# 
Solutions obtained for values of sn greater than one require that the 
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colunn pass through a ro 'ion of instab t lit and aro thus possi le only 
in thojr ' or wit' the use of external constraints.) 



Figure 3. Schematic Represent at ion of the Experimental Set-up 

Figure 3 is a schematic representation o: 1 the experimental set-up. 

If the piston is frictionless then P may be replaced by its equivalent 
in terms of internal pressure jfa * The fluid can not sustain shear 
so that the modulus of elasticity and the second moment of cross-sectional 




B 



area depend only on the material and cross-section of the cylinder. 
Plaoing the equation in terms of the internal pressure we obtain 

7t 1 E At 

( i)‘A, ' 

Where A Q is the area of the cross-section, and C is the radius of 
gyration of the oross-section of the cylinder# 
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Ill Description of Apparatus 
1. Selection of Material 

The buckling load, and thus the internal pressure producing buckling, 
is directly proportional to the modulus of elasticity# Aluminum alloy 
has the lowest modulus of elasticity among the readily available tubing 
materials and was therefore choson# 

2# Choice of Geometrical Properties 

The principal limitation placed on the geometric properties was that 
the tube buckle before the hoop tension reached the yield point. A 
slenderness-ratio of 100, which fulfills this requirement, was selected 
as the minimum to be tested. 

Establishing a minimum slenderness-ratio and a maximum internal 
pressure (see Section 5) dictates a maximum diameter-to-thickness ratio. 
Among the geometrically similar tubes which oould be used outer diameters 
of 3/4, 1 and l-l/4 inch were selected because they provide lengths whioh 
could be conveniently handled. A thickness of .035 inch was required for 
the l-l/4 inch tube. This thickness also provided satisfactory charac- 
teristics for the other two sizes. 

3. End Fittings 

Fig. 4 is a schematic diagram of the end fittings. The tube and 
end fittings are joined by extra fine threads, 56 threads per inch. The 
upper end fitting has a cylinder joined to it by threads. (Henceforth 
this will be considered one piece and called n the upper fitting.”) The 
upper portion of the piston, which contacts the horizontal surface of 

the testing machine, is a spherical surface (radius 2©75 inches). The 

» 

center of this surface is on the axis of the piston in the plane of the 
face of the upper fitting. (Face of a fitting is defined as that flat 
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s t’ rou fh w : ch the tube is t* roadod. ) 





b) Upper 



Fi ;ure 4. bnd Fittings 

The lower fitting is similarly formed with the center of the 
spherical surface (radius 1*75 inches) on its axis in the plane of its 
face. The lower fitting has a passage, fitted with a 1/4 inch pipe 
nipple, to allow entry of the hydraulic fluid. 

The effect of the spherical surfaces is to cause the column to act 
as if it were in the hinged-end conf iguration with an effective length 
equal to the distance between the faces© 

The piston and upp r fittings (Fir, 5-a) retain the hydraulic 
fluid and dictate the axial load transmitted to the tube 0 The piston 
is fitted with an O-rinr seal to limit leakage. The same piston is used 
in testing the threo different sizes of columns. This results in a 
positive load-ratio, a zero load-ratio or a negative load-ratio; 
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a) Piston, and upper 
fitting with 1 
inch tube 




b) Upper fitting 
with 3/4 inch 
tube 




c) Upper fitting 
with 1-1/4 inch 
tube 



Fi uro 



C'.’lind'T, Pist'ir, 



and Fni Fitting 
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depending on whether the inner diameter of the tube is less than 
(Fig, 5-b), equal to (Fig, 5-a), or greater than (Fig, 5-o ) the diameter 
of the piston, (The diameter of the piston is taken to be the inner 
diameter of the cylinder,) 

4, Deflection Measurement 

Lateral deflections were measured by means of two dial gages set at 
right angles to each other and mounted, at the mid-point of the column, 
in a plane normal to the axis of the undeformed column, 

5, Pressure Application 

The pressure was applied by a dead-weight gage testing machine, 
using hydraulic fluid. The apparatus was designed to carry pressures to 
2000 psi. It was felt that pressures greater than this would complicate 
the sealing of the threaded section and require a more complex 0-ring 
design. 

The fluid is oarried from the gage tester to the lower fitting 
through 10 feet of hose having an inner diameter of 3/l 6 inch and an 
outer diameter of 33/64 inch. The hose consists of a synthetic rubber 
inner tube, a single wire braid reinforcement, and a synthetio rubber 
cover. 

Fig, 6 shows a photograph of the apparatus assembled and ready for 

test. 
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Figure 6. Photograph of the Apparatus 
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IV Experimental Procedure 

In deriving the theory it was assumed that the piston was friction- 
less# The experimental procedure was devised to eliminate the effects 
of friction between the piston and cylinder# The first method attempted 
was to apply the pressure, then turn the piston within the cylinder# 

This method appeared to eliminate the friotion, but the resulting lateral 
movement of the column made the lateral deflection readings unreliable. 

The method finally used to eliminate the frictional effect was to apply 
the pressure and then adjust the testing machine to a load equal to the 
product of the piston area by the applied pressure* The column was 
aligned in a vertical position by using the plumb bubbles on a 28-inoh 
oarpenter’s level# 

The sequence of operation was as follows: 

1) The column was plaoed in a vertical position by 
means of the level; 

2) a fraction of the oritioal pressure was applied 
by loading the gage tester; 

3) the load on the testing machine was adjusted to 
that load which it would oarry if the piston were 
friotionless; 

4} the lateral deflections were read from the dial gages; 

5) the piston displacement was oheoked to insure that it 
had not reached the end of its travel; 

6) the pressure applied by the gage tester was increased# 

Steps 3) through 6) were repeated until the deflection read in Step 4) 

/ 

reaohed the limit of the dial gage. The oritioal pressure was 
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obtained usin^ a nothod su^osted by Southwell (/Vppendix I 



f Section l). 
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V Results and Conclusions 



1. Tabular Results 

The theoretical and experimental values of oritioal pressure and 
unit load (unit load is the quotient of the total axial load on the 
column by the cross-sectional area of the cylinder) are summarised in 
Table 1* Additional data appears in Appendix I, Section 4# 

2# Graphical Results 

The theoretical curve of unit load versus slenderness-ratio, with 
experimental points plotted, is given in Fig# 7# 

3# Conclusions 

The experimental data indicate that the equation of Euler 

7T* E I 

V 

is applicable to hollow cylindrical shapes subjected to internal 
pressure, where: P is the total axial load on the oolumn; and the 

modulus of elasticity and second moment of cross-sectional area depend 
solely on the material and cross-section of the cylinder# 
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Table 1 



Column 


Slenderness 


Critical 


Pressure (psi) 


Unit 


Load 


% Dif- 


Number 


ratio 


Theoretical Experimental 


Theoretical Experimental 


ference 






3/4 


inoh tube 








1 


105.1 


1092 


1074 


9460 


9300 


-1.7 


2 


105.1 


1093 


1091 


9460 


9450 


-0.2 


3 


120.2 


835 


820 


7220 


7090 


-1.9 


4 


120.2 . 


835 


819 


7220 


7080 


-1.9 


5 


147.8 


555 


543 


4790 


4690 


-2.2 


6 


147.8 


553 


543 


4790 


4700 


-1.8 


7 


202.6 


296 


296 


2550 


2550 


0 


8 


202.6 


297 


292 


2550 


2510 


-1.7 






1 


inoh tube 








9 


100.5 


1643 


1625 


10360 


10250 


-1.1 


10 


100.5 


1646 


1649 


10360 


10380 


/0.2 


11 


118.7 


1172 


1189 


7430 


7540 


/1.4 


12 


118.1 


1190 


1194 


7500 


7520 


■/0.3 


13 


142.2 


816 


800 


5180 


5070 


-2.0 


14 


142.3 


816 


811 


5170 


5140 


-0.6 


15 


185.6 


474 


465 


3040 


2980 


-1 0 9 


16 


185.5 


475 


474 


3040 


3030 


-0 o 5 






1-1/4 


inoh tube 








17 


. 101.5 


2068 


2012 


10160 


9890 


-2.7 


18 


101.5 


2069 


2024 


10160 


9940 


-2.2 


20 


118.0 


1508 


1503 


7510 


7480 


-0.3 


21 


136.3 


1138 


1120 


5630 


5540 


-1.6 


22 


136.3 


1136 


1118 


5630 


5550 


-1.6 


23 


152.4 


917 


917 


4510 


4510 


0 


24 


152.4 


919 


917 


4510 


4510 


-0*2 



18 




e 



• ?'■ errot eft! "'irre w * . Jxpor ir.O'. 1 nl I’o^ts 



19 



BIBLIOGRAPHY 



1. L. Euler, On the Strength of Columns, as translated by J • A. 
Van Den Broek, American Journal of Physics, Vol. 15, No# 4, 
309-318, July-August 1947# 

2. R# V. Southwell, On the Analysis of Experimental Observations 
in Problems of Elastic Stabilit ', Proceeding of the Royal 
Society, London, Series A, Vol. 135, P. 601, 1932. 

3# S# Timoshenki, Theory of Elastic Stability, McGraw-Hill Book 
Canpany, Inc., 1936. 



20 



APPENDIX I 



« 

ANALYSIS OF EXPERIMENTAL DATA 



1 • Southwell's Method 

In the ideal case of column buckling there would be no lateral 
deflection up to the critioal value of the load. At the critical load 
the column would then begin to deflect. Owing to various kinds of imper- 
fections, such as some unavoidable initial curvature of the column, 
eccentricity in application of the load and non-homogeneity of the 
material, the column begins to deflect with the beginning of loading and 
usually fails before Euler's load is reaohed. 

A useful method of determining the critioal load from the test data 
was proposed by R. V. Southwell [2 J • Assume that the deflection of a 
oolumn, under a load that is below the critioal value, is due to initial 
curvature and eccentricity of load application. If the deflection is 
expressed as a trigonometric series, when the load approaches the critical 
value the first term in the series predominates. 

If the deflection d is represented by the first torn of the series, 
the following linear relation exists 

d- & i- -a 



where is a constant. A cartesian plot of d versus yields a curve 

r 

whose slope asymptotically approaches the critical load. 

The experimental data of Column 11 are plotted in Figs. 0 through 
12. The nature of the curve and its asymptotic behavior can be observed# 
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2o Sample Calculations 

The oritioal load used in the theoretical derivation is tho product 
of the internal pressure and the piston area* In the experimental 
prooedure it was advantageous to consider the internal pressure to be 
the principal variable* For this reason the data are reduced in terns of 
the internal pressure* 
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Column 11 


Run 1 










Gage 1 




Gage 


2 




4 


d 


d /-fc 


d 




d //. 


200 


2.5 




5 




.025 


400 


3.5 




13 




.032 


600 


4.o 




22 




.037 


700 


5 




32 




.046 


800 


6 




42 




.054 


850 


6 




49 




.058 


900 


6.5 




63 




.070 


950 


6.5 




72 




.076 


975 


6 




84 




.086 


1000 


6 




100 




.100 


1025 


6 




108 




.105 


1050 


4 




133 




.127 


1075 


4.5 




160 




.149 


1100 


0.5 




187 




.170 


1125 


-8 




250 




.222 


Gage 1 shows a 


possible variation of the 


data. 


At moderate 



pressures sene term of the series other than the first predominates. 
As the critical pressure is approached the principal mode of 
distortion predominates. The data are not plotted since no trend 
can be determined# 

Data from Gage 2 are plotted in Fig # 8. 
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Column 


11 


Run 2 










Gage 


1 








Gage 2 




* 


d 






<vk 




d 




200 


1 






.0050 




3 


.015 


400 


3.5 






0 OO 88 




8 


.020 


600 


7.5 






.0125 




18 


.030 


700 


9.5 






.0136 




24 


.034 


800 


15 






.0187 




34 


.042 


850 


16.5 






.0194 




41 


.048 


900 


18.5 






.0210 




49 


.056 


950 


22.5 






.0237 




64 


.067 

« 


975 


23.5 






.0241 




69 


.071 


1000 


26.5 






.0260 




83 


.083 


1025 


29.5 






.0288 




96 


.093 


1050 


32.5 






.0310 




110 


.105 


1075 


36.5 






.0340 




129 


.120 


1100 


39.5 






.0360 




155 


.141 


1125 


41.5 






O 0369 




173 


.154 


1150 


47.5 






.0413 




250 


.218 


Data 


from Gage 


1 


ere 


plotted 


in 


Fig. 9. 




Data 


from Gage 


2 


are 


plotted 


in 


Fig. 10. 





25 



Column 11 Run 3 



A 


Gage 1 
d 


d/yb 


Gage 2 
d 


d /V> 


200 


0 




5 


.025 


400 


1 




12.5 


.031 


600 


1 




23.5 


.039 


700 


0 




32.5 


.046 


800 


0 




43.5 


.054 


850 


0 




53.5 


.063 


900 


0 




63.5 


.071 


950 


-1 


.0015 


73.5 


.078 


975 


-3 


.0030 


88.5 


.091 


1000 


-3 


.0030 


93.5 


.094 


1025 


-5 


.0049 


106.6 


.104 


1050 


-8 


.0076 


138.5 


.132 


1075 


-11 


.0102 


164.5 


.153 


1100 


-15.6 


.0141 


171.5 


.156 


1125 


-26 


.0231 


236.5 


.210 


Data from 


Gage 1 are 


plotted 


in Fig. 11. 





Here is an example of the first term of the series 
predominating as the critical pressure is approached. 
In this oase a trend is established. 

Data from Gage 2 are plotted in Fig. 12. 
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3. Criteria Used in Evaluating Critical Pressures. 

Let the points obtained from the data be numbered "1" through "n" 
in order of inoreasing load. For a "good" data set the following 
generalisation oan be made: 

1) The slope of tho line joining points "n-2" and "n" will 
be approximately equal to the slope of the lino joining 
points "n - 1" and "n." 

2) The best estimate of the asymptotic slope, without 
extrapolation, is the 3lope of the line joining points 

"n - 1" and "n." 

Table 2 lists the slopes obtained from the last three points of 
Figs. 8 through 12 

Table 2 

Slopes Obtained from Data Sets of Column 11 

Run 1 



Gage 1 




Gage 


No ourve 




1231 






1210 




Run 2 




Gage 1 




Gage 


1510 




1233 


1362 


Run 3 


1202 


Gage 1 




Gage 


1163 




1262 


1165 




1202 
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The following criteria were set up to evaluate the information 
obtained frcm the plots, Data sets in which the deflection ohangea 
dirootion within the last three readings were discarded without plotting. 
Data sets in whioh the two slopes oanputed from the last three points 
differed by more than five per cent wore discarded. Data sets in whioh 
the two slopes agreed within five per oent were assigned a weight of 
unity. Data sets in whioh the two slopes were within one per oent were 
assigned a weight of 2, The weighted average of the slopes of the lines 
joining the last two points was taken as the critical pressure of the 
column. 

The above system of evaluation, though it does not prevent discard- 
ing good information, makes it unlikely that poor information will be 
included, 

4, Tabulation of geometrical properties and oritioal pressures. 

The geometric properties, oritioal pressures, aid average critical 
pressure for each column are given in Table 3, Method of determining 
theoretical values used in Table 3 may be found in the section of 
Appendix I indicated: 

Geometric Properties Seotion 5 

Critioal Pressures Section 3 

Average Critical Pressure Section 3 
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Table 3. Tabulation of Geometric Properties, Critical Pressures 
and Average Critical Pressure for Runs Used 
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5. Measured Quantities 

1 • The length of the column was measured with a 6 foot steel tape. 

It is expected that measurements are within *03 inch. 

2. The outer diameter and thicVmess were measured with a micrometer. 
It is expected that measurements are within .0005 inch. 

3. The radius of gyration was obtainod analytically using the outer 
diameter and thiokness. 

4 0 The weight of the column was obtained to the nearest 0.1 ^ran. 

5. The area was obtained analytically by dividing the tube weight 
by the product of the length by the density# The tabulated density of 
24S-T3 aluminum alloy is 0.100 pounds per oubio inch. 

6# The tabulated value of the modulus of elasticity of 24S-T3 
aluminum alloy is 10#6 x 10^ lbs per inch* This value was verified 
exper iment al ly • 

7. The theoretical critical pressure was computed using a calcu- 
lator. An excess number of figures was carried to avoid round-off 

errors# 
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